Making the best of mixed-field orientation of polar molecules: A recipe for achieving 
adiabatic dynamics in an electrostatic field combined with laser pulses 
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We have experimentally and theoretically investigated the mixed-field orientation of rotational- 
state-selected OCS molecules and we achieve strong degrees of alignment and orientation. The 
applied moderately intense nanosecond laser pulses are long enough to adiabatically align molecules. 
However, in combination with a weak dc electric field, the same laser pulses result in nonadiabatic 
dynamics in the mixed-field orientation. These observations are fully explained by calculations 
employing, both, adiabatic and non-adiabatic (time-dependent) models. 



PACS numbers: 37.10.Vz, 33.15.-e, 33.80.-b, 42.50.H2; 

Creating oriented samples of polar molecules, i.e., 
molecules with their dipole moment preferentially point- 
ing towards one hemisphere rather than the opposite, 
has been a long standing goal in molecular sciences. It 
was originally motivated by the crucial role played by 
orientation in chemical reaction dynamics [T]. More re- 
cently, its importance in novel applications such as (fs- 
time-resolved) photoelectron angular distributions 
diffraction- from- within [5], or high-order harmonic gen- 
eration [6] has been realized. 

Early methods exploited purely electrostatic fields. Us- 
ing an electric multipole focuser, molecules in a single 
low-field-seeking quantum state can be selected due to 
their first-order Stark effect The degree of orien- 

tation is determined, and also limited, by the selected 
state. Alternatively, a strong homogeneous electric field 
can create so-called brute- force orientation [lOl E] . This 
method requires very high electric field strengths and 
works best for rotationally cold molecules with large per- 
manent dipole moments. 

In 1999 a method based on the combined action of a 
moderately intense, nonresonant laser field and an elec- 
trostatic field was proposed [T2]. For the case that the 
laser field is turned on significantly more slowly than the 
rotational period(s) of the molecule adiabatic behavior 
was assumed. The time-independent calculations showed 
that the degree of orientation could be nearly perfect 
under conditions present in many experimental setups. 
Furthermore, the degree of alignment, i.e., the confine- 
ment of the molecular axes to space-fixed axes, could 
also be very high. In addition, the method should be 
generally applicable to a broad range of molecules and, 
therefore, promises the availability of strongly oriented 



and aligned molecules for various applications. Experi- 
ments performed in the first half of the 2000s showed the 
feasibility of the method but the degree of orientation 
observed was moderate [131 [14]. A major reason for the 
weak orientation was that while the individual pendular 
states are strongly oriented, these states arise in pairs 
whose members are oriented oppositely with respect to 
one another. Consequently, the resulting overall degree 
of orientation, obtained as the weighted average over the 
populated quantum states, diminishes compared to what 
is expected for very cold or even single-state molecular 
ensembles. A significant improvement in the experimen- 
tal capabilities was reported in 2009 when quantum-state 
selected molecules were employed as targets leading to 
strongly enhanced orientation [TSHiT] . However, it was 
already realized that an adiabatic description is not suf- 
ficient to reproduce the experimental observations [18j . 

In the present work, we seek the maximum of achiev- 
able orientation, as predicted by the original adia- 
batic description. Therefore, we prepare a nearly 
pure (92^5 %) rotational-ground-state ensemble of OCS 
molecules [12] and use a laser pulse that is sufficiently 
strong to ensure sharp alignment and that is turned 
on a 100 times slower than the rotational period of 
the molecules. Our experimental observations are, how- 
ever, at odds, both qualitatively and quantitatively, 
with the predictions of the original theory [T^]. In- 
stead, the experimental findings, exploring the depen- 
dence of the orientation on both the laser intensity and 
on the static field strength, can be rationalized by solving 
the time-dependent Schrodinger equation describing the 
mixed-field orientation. Our analysis directly shows the 
nonadiabatic coupling of the two sublevels of the near- 
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FIG. 1. (Color online) (a) Schematic of the field configura- 
tions showing the polarization directions of the alignment and 
the probe pulses as well as the static field inside the VMI spec- 
trometer and the definition of the angle /3. (b),(c) S"*" ion im- 
ages for /3 = 45° and 135°, respectively and Egtat = 571 V/cm. 
The rings in the image depict the limits for S"*" ions from the 
S''' + CO^ channel used in the calculation of the degree of 
orientation. 



degenerate doublets, created by the laser field, and allows 
us to predict the experimental conditions needed to en- 
sure adiabatic dynamics. 

The experimental setup has been described in detail 
before [T71 [5^ and only a few important details will 
be pointed out, see Fig. J| a). A pulsed molecular beam 



is formed by expanding a mixture of 1 mbar of OCS and 
10 bar of neon into vacuum through a pulsed valve. The 
molecular beam is skimmed twice before entering a 15- 
cm-long electrostatic deflector. Here it is spatially dis- 
persed in the vertical direction according to the effec- 
tive dipolc moments of the quantum states [T7]. Here- 
after, the molecules travel into a velocity map imag- 
ing (VMI) spectrometer where they are crossed by two 
pulsed laser beams. The first pulse (Eaiign, A = 1064 nm, 
tfwhm = 8 ns, linearly polarized) provides the laser field 
for the mixed-field orientation whereas the weak static 
field, Estat, exploited for the orientation is (inherently) 
provided as part of the VMI spectrometer, which also de- 
fines its direction. The second pulse (probe, A = 800 nm, 
tfwhm — 30 fs, linearly polarized) is used to charac- 
terize the orientation and alignment by multiply ioniz- 
ing the molecules, this is followed by Coulomb explo- 
sion and imaging of the recoiling S+ fragments on a two- 
dimensional detector. 

The strongest orientation is expected when Eaiign is 
parallel to Egtat- This geometry is, however, not well 
suited for the ion imaging method to characterize the 
orientation because the experimental observable, the S''" 
ions, will then be localized in the center of the detec- 
tor. Consequently, all measurements are conducted with 



E; 



align_ 



rotated by an angle (3^0 with respect to Egtat 



(see Fig. 1 



a)). Figures[TJb) and (c) show examples of S+ 
ion images for /? = 45° and 135° (equal to -45°). The S+ 
ions from the Coulomb explosion channel S+ -I- CO^, ap- 
pearing in the outermost part of the images, are highly di- 
rectional and provide direct information about the align- 
ment and orientation of the OCS molecules at the time 
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FIG. 2. (Color online) Orientation ratio for j3 = 45° as a func- 
tion of laiign, for the weak (a) and the strong (b) static field, 
showing the experimental results (black circles), the adiabatic 
calculations (green dotted line) and the time-dependent cal- 
culations (blue solid line). 



of ionization. 

The strong angular confinement of the S+ ions shows 
that the OCS molecules are sharply one-dimensionally 
aligned along Eaiign- In addition, a pronounced asym- 
metry of the S"*" ions emitted either along or opposite 
to Egtat, with an excess of S"*" in the upper (lower) re- 
gion for /3 = 45° (135°) is observed. This shows that 
the molecules are oriented with the S-end preferentially 
pointing toward the detector screen - as expected and in 
agreement with previous studies [2T] . To quantify the de- 
gree of orientation only ions from the S"*" -I- CO"*" channel 
are considered. We then specify the orientation by the 
ratio Nup/Ntot of the number of these ions in the upper 
half of the image Nup compared to the total number of 
ions Ntot, from this channel. 

In the first set of measurements the degree of orien- 
tation is recorded as a function of the alignment pulse 
intensity, laiign, for two values of Egtat- The results are 
For low values of laiign the orienta- 



shown in Fig. 2 



tion ratio is almost the same for the two static fields but 
for laiign > 2.5 X 10^^ W/cm^ the results differ. For the 
strong static field the orientation reaches a maximum of 
approximately 0.8 at laiign = 5 x 10^^ W/cm^ and re- 
mains essentially constant out to 1.4 x 10^^ W/cm^. In 
contrast, for the smaller static field the maximum ori- 
entation occurs already at laiign = 3 x 10^^ W/cm^ and 
the degree of orientation decreases as laiign is further in- 
creased, dropping to 0.70 at laiign = 1-4 x 10^^ W/cm^. 
The calculated degree of mixed-field orientation using the 
adiabatic model for /3 = 45°, as in the experiment, is 
shown in Fig. 2 Here, we have used rotational-state 



populations, in the coordinate system of the electric field 
in the deflector, with 92 % in the |0,0) state, adiabat- 
ically corresponding to the field-free J = M = state, 
4% in the state and 4% in the state [12]- 

These states are projected onto a coordinate system for 
the mixed-field orientation that is defined by Eaiign- The 
properly symmetrized states are then |0,0, e), |l,l,e), 
and |1, l,o), where e and o denote even and odd parity 
with respect to the plane defined by Egtat and Eaiign, 
respectively. The volume effect ITS] is accounted for by 



3 




Experimental 
.E,tat =286V/cm 
■ Estat =571V/cm 



Time dependent 

Estat =286V/cm 

Estat =571V/cm 



Adiabatic 

Estat =571V/cm 



FIG. 3. (Color online) The orientation ratio as a function 
of P for laiign = 9.1 X lO'^^W/cm^ for the experiment using 
Estat = 286V/cm (black circles) and Estat ~ 571V/cm (red 
squares), for the adiabatic calculations, which are identical 
for the two static fields (green dotted line), and for the time- 
dependent calculations using Estat ~ 286V/cm (blue dashed 
line) and Estat = 571 V/cm (blue solid line). 



using an experimentally determined cubic dependence on 
the probe pulse intensity. These calculations predict that 
the orientation is independent of the applied static field. 
Following a rapid initial rise it reaches a value of 0.96 al- 
ready at 4 X 10^^ W/cm^ and remains constant. Clearly, 
these predictions are at odds with the experimental find- 
ings: The simulated degree of orientation is much too 
strong and it does not reproduce the decrease of the ori- 
entation at intensities above 3 x 10^^ W/cm^, that is ex- 
perimentally observed for the case of the smaller static 
field. 

the 



In the second set of measurements, shown in Fig. 3 



degree of orientation is recorded as a function of /3 for the 
weak and the strong static fields and for a fixed value of 
laiign of 9.1 X 10^^ W/cm^. For both static field strengths 
Nup/Ntot decreases monotonically as /3 increases from 30° 
to 150°. At all f3 values the strong field leads to stronger 
orientation than the weak field. The N^p/Ntot ratios cal- 
culated from the adiabatic model are essentially identical 
for the two static field strengths. The sharp rise (fall) of 
the curve to a value close to 0.96 (0.04) as /? is increased 
(decreased) below (above) 90° shows that very strong ori- 
entation is reached already for a very modest static elec- 
tric field along Eaiign- This calculated behavior of the 
orientation differs qualitatively as well as quantitatively 
from the experimental results. 

To obtain a better model of the mixed-field orienta- 
tion process we solve the time-dependent Schrodinger 
equation using the experimental field configurations 
and rotational-state populations [22]. The results for 
Nup/Ntot as a function of laiign are shown in Fig. 2 



The predictions of stronger orientation for the strong 
static field and, in the weak static field case, the de- 
creasing orientation at increasing intensity for laiign > 
3 X 10^^ W/cm^ are in line with the experimental find- 
ings. Moreover, the smooth /3-dependence of the orien- 
tation, shown in |Fig. 3[ is fully captured by the time- 
dependent calculations. Quantitatively, the calculated 




- |0,0,e) 
-|l,l,e 

■'|l,0,e>p_ 
-|2.2,e>p 



2 3 4 5 6 
Ialign(10" W/cm^) 
t(ns) 





2 3 4 5 6 

Ialign(10" W/cm^) 

FIG. 4. (Color online) (a) Energy of the 4 lowest lying rota- 
tional eigenstates as a function laiign (time). The inset shows 
the relevant energy and intensity ranges where the formation 
of the near-degenerate doublets occurs, (b) Polar plot repre- 
sentation of the wavefunctions for the 4 states shown in (a) at 
laiign = 9 X 10" W/cm^. The single-headed arrow shows the 
direction of the static field, (c) The squares of the coefficients 
for the projection of the time-dependent pendular wavefunc- 
tion of the absolute ground state onto the adiabatic pendular 
state basis (|0,0,e) (black) and |l,l,e) (green)) as a func- 
tion of laiign (time) for /3 = 30° (solid line), 45° (dotted line) 
and 89.5° (dashed line). Estat = 286 V/cm for all data. 



values overestimate the degree of orientation. This could 
partly be due to temporal substructure in the experi- 
mentally applied laser pulses, which could induce more 
nonadiabatic population transfer. Overall, these non- 
adiabatic-model calculations are in much better agree- 
ment with the experimental results than the adiabatically 
calculated ones. 

The underlying physical picture for understanding the 
failure of the adiabatic model is obtained by considering 
the evolution of the states during the turn-on of the align- 
ment pulse. Before the pulse, the rotational states are es- 
sentially described by field-free rotor states. As the laser 
field strength increases, the states are hybridized by the 
combined action of the laser and static fields. In Fig. 4[ a) 
the formation of doublets of nearly-degenerate pendular 
states in the strong laser field regime is shown. For the 
laser and static fields used in the experiment the abso- 
lute ground state |0, 0, e)p is right-way oriented, i. e., the 
permanent dipole moment is pointing along Egtat- The 
upper level, | 1, 1, e)^ , of the lowest doublet is wrong- way 
oriented - see Fig. 4lb). 



As the alignment field is turned on the states that even- 
tually form the near-degenerate doublet are coming closer 
together. This is illustrated in Fig. 4[ a) for the |0,0, e), 
|l,l,e) pair and the |1,0, e), |2, 2,e) pair and it results 
from the ac Stark interaction. When the energy splitting 
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the |0,0, e)p and |l,l,e)p states 
which is the case for a weak static electric field. 



within a pair approaches the coupling strength due to 
the dc Stark interaction between the two sublevels, the 
two states in each pair can mix because they have the 
same symmetry, provided the laser and static fields are 
non-perpendicular. This will result in population trans- 
fer between the oriented and anti-oriented states. The 
probability for mixing, corresponding to a crossing from 
one state of the doublet to the other, is determined by 
the rate of the turn-on and the energy separation between 
1, 1, e)„ states. If the splitting is small, 

the two 

states will be strongly coupled. To ensure fully adiabatic 
transfer it is necessary to turn on the laser field on a 
timescale slower than the inverse of the energy splitting 
of the near-degenerate doublet. This time can be much 
longer than the rotational period of the molecule. For the 
lowest doublet formed in OCS with Egtat — 286 V/cm and 
laiign = 9.1 X 10^^ W/cm^ our calculations show that the 
alignment pulse must be 50 ns long to ensure adiabatic 
transfer. The 8 ns pulses used in the current experiments 
do not fulfill this adiabaticity criterion although they are 
a hundred times longer than the rotational period of OCS 
- the condition previously considered sufficient for adia- 
batic behavior. 

The population transfer is illustrated in Fig. 4| c), 



where the decomposition of the time-dependent state, 
which starts as the rotational ground state, in terms of 
the pendular states, is shown during the time interval 
representing the turn-on of the alignment pulse for three 
/3 values. The field-free ground state is not transferred 
adiabatically to the pendular ground state: For /3 — 45°, 
the final state is decomposed into 74.06% |0,0, e)p and 
25.94% |1, l,e)p. Therefore, the resulting degree of ori- 
entation falls below that expected for a pure adiabatic 
transfer since the |l,l,e)p state is wrong-way oriented. 
Similarly, other field-free rotational states are mixed with 
different pendular states during the turn-on, for instance 
the initial state ends up in a superposition of 

13.00% |0,0,e)p, 37.10% |l,l,e)p, 35.64% |l,0,e)p, and 
14.26% |2,2,e)p. For (5 = 89.5°, the electric field along 
the molecular axis is small, and the |0, 0, e) and |1, 1, e)p 
states contribute with 50.32% and 49.68%, respectively, 
to the time evolution of |0,0), resulting in a vanishing 
orientation. By contrast, alignment is expected to re- 
main strong since both the |0, 0, e)p and |l,l,e)p states 
imply tight confinement of the molecular axis along the 
laser field polarization, see Fig. 4[ b). The experimental 
observations for perpendicular fields do indeed show no 
orientation but strong alignment [23j . 

In summary, the combined action of a moderately 
strong laser field and a weak electrostatic field re- 
mains an attractive approach for creating tightly oriented 
molecules, but to fully exploit the potential of the method 
it is necessary to redefine the meaning of adiabatic condi- 
tions. Unlike alignment, where adiabaticity is ensured by 
turning on the laser field slower than the rotational pe- 



riod of the pertinent molecule, adiabatic transfer in orien- 
tation necessitates that the laser field be turned on slower 
than the inverse of the minimum spacing between the two 
pendular states in a doublet. This has repercussions for 
designing experimental parameters such that the degree 
of orientation be optimized. In the case of OCS, our cal- 
culations show that when Egtat — 286 V/cm and f3 = 45° 
adiabatic transfer of the |0, 0) state to the |0, 0, e)p state 
is obtained if a transform-limited laser pulse with a Gaus- 
sian pulse duration (full width half maximum) of 50 ns 
is used. In practice, such pulses are not easily supplied 
by lasers typically present in laboratories. Considering 
instead the 10 ns output from the widespread Nd:YAG 
lasers adiabatic transfer will occur for a static electric 
field of 2kV/cm: The increased static field leads to a 
larger minimum spacing of the doublet and, thus, relaxes 
the requirement for the slowness of the laser-field turn-on. 
Such static fields are compatible with, for instance, VMI 
spectrometers. For a pulse durations of 500 ps, which 
would be relevant for using the stretched output from 
amplified Ti-Sapphire lasers, a static field of lOkV/cm 
is needed to ensure adiabatic conditions. Again this is 
compatible with electron or ion spectrometers |TB) . 

We note that the lack of adiabatic behavior will also be 
influenced by avoided crossings between rotationally ex- 
cited states at low laser intensities. For larger molecules, 
where the rotational level structure is quite complex, this 
effect is expected to be particularly important [T51 HI] , 
but an increase of the static field strength should enhance 
the degree of orientation as already demonstrated experi- 
mentally for several asymmetric top rotors [2 [13 |T71 1^ . 
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